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oa Shaft CenterLINES 


Estimating vibration at critical 


locations using available measurements 


and machine configuration data 


Nomenclature 

AS = Element shaft cross sectional 
area 

Ad = Disk cross sectional area 

Db = Bearing damping 

li = Element shaft inside diameter 


bending moment of inertia 

lo = Element shaft outside diameter 
bending moment of inertia 

Kb = Bearing stiffness 

Kd = Bearing direct dynamic stiffness 

Kq = Bearing quadrature dynamic 
stiffness 

Ld = Disk length 

Ls = Element shaft length 

SNR = Signal to noise ratio 

A = Average circumferential fluid 
velocity 

pd = Disk density 

ps = Element shaft density 


tl = Forcing function rotative speed 
2 = Rotor rotative speed 


ne of the problems associated 

with establishing acceptable 

vibration levels for rotating 
machinery is noncollocation of critical 
vibration and available measurement 
locations. Noncollocation of measure: 
ment and critical machinery monitoring 
locations has been a problem faced by 
machinery engineers since the first mea- 
surement probes were installed on 
machinery. Critical machinery monitor- 
ing locations are any location along the 
rotor with a reduced clearance or a high 
potential for rotor to stator contact. 

It is difficult to estimate what the mea- 
sured vibration at the bearings will be 
just before the available clearance is 
exceeded at an internal location. 


Unknown forcing functions, operation 
through several lateral modes, and scar- 
city of measurement locations compli- 
cate the estimation process. For rigid 
machines, these problems are easily 
overcome. However, for flexible 
machines, which must pass through sev- 
eral lateral modes getting to operating 
speed, this becomes more difficult. The 
difficulty increases because the correla- 
tion between measured and critical 
vibration responses changes with imbal- 
ance distribution, operating speed, nat- 
ural rotor modes, and other operating 
conditions. The creation of two mea- 
surement locations near the bearings, 
namely, “modal” probes, instead of just 
one, increases the measurement loca- 
tion density, which can lead to an 
improvement in the estimation of the 
actual operating mode. If a multi- 
hearing machine train (Figure 1) is con- 
sidered, several data interpolation/ 
extrapolation techniques could be used, 
some of which are outlined below. 

Linear Extrapolation: For linear 
extrapolation, the bearing closest to the 
desired location is selected. The equa- 
tion of the straight line through the two 
data points on either side of the bearing 
can then be used to find the response 
amplitude and phase at other axial loca- 
tions. This technique is useful for rigid 
modes, but starts producing errors as 
soon as the rotor begins to bend. 
Clearly, another technique must be 
found if extrapolation for rotor bending 
modes is required, 

Polynomial Extrapolation: Poly- 
nomial extrapolation requires that a 
polynomial of order one less than the 
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number of measurement locations be 
generated. This polynomial equation 
can then be used to calculate the vibra- 
tion response at any axial location. This 
technique might seem to be ideal for the 
task at hand. However, high order inter- 
polation polynomials tend to produce 
large errors for a great many functions, 
especially if the measured signals are 
contaminated by noise, 

Spline Fits: The mathematician's 
answer to the errors introduced by high 
arder polynomials are spline fits. Gener- 
ally, spline fits replace a single high 
order polynomial with several lower 
order polynomials, each of which is used 
for some portion of the data range. To 
smooth the transitions from one poly- 
nomial to the next, both the response 
amplitude and phase, and several of its 
derivatives, must be identical to adjoin- 
Ing polynomials at the transition points. 
The measurements give the response 
vectors, but there are no measurement 
devices that measure derivatives 
directly. Therefore, they must be calcu- 
lated from the response amplitudes. 
Modal probes allow the first derivative to 
be approximated reasonably well. How- 
ever, to calculate the second derivative, 
another point is necessary. A spline fit 
algorithm would go to the next available 
data point to get this required informa- 
tion. However, if we notice how far that 
point is from the desired location, across 
a rotor span, you have to wonder 
whether this point provides any useful 
information about the seconecl derivative 
at the original location. Because of this, 
spline firs are not good candidates for 
extrapolation of responses along 
machine trains. 
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Cubic Polynomials: If cach rotor 
span is processed separately, a cubic 
polynomial could be constructed which 
passes through the four measurement 
locations at the support bearings. This 
allows good approximations up to and 
including the second bending mode for 
the selected rotor span. Assuming that 
other, more relevant, information is not 
available, this technique probably gives 
the best results in the field. 

However, some additional informa- 
tian exists which could lead to better 
extrapolated responses. The physical 
configuration of the system is useful data 
that is often overlooked, This includes 
the bearing locations and their stiff- 
nesses, the operating speed and the 
physical description of the rotor. There 
are several methods suitable for using 
system information and measured data 
to extrapolate vibration responses. 
Eigenvectors can be computed and used 
as anew modal coordinate system. The 
number of degrees of freedom can then 
be reduced by discarding all coordinates 
except those that describe the desired 
modes. Either optimization or autocor- 
relation techniques can be used to deter- 
mine the coefficients of each mode. The 
response at the desired location can 
then be obtained by translating back to 
the original coordinate system. 

An alternate way to compute vibration 
responses at the measured locations is 
lo optimize the imbalance force distribu- 
tion to minimize the error between 
these computed responses and the 
actual measured responses. This alter: 
nate approach is based on finite element 
or transfer matrix techniques, Similar 
techniques have been used for machin- 
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ery parameter identification by other 
researchers (Nordmann, Diewald 
(1990), Mottershead (1991), Muszynska, 
et al (1992)). A computer program 
implementing this process, along with 
preliminary results from experimental 
testing, is presented. 


Finite element optimization 
technique 

The flow diagram for the computer 
program is shown in Figure 2. The 
machinery engineer must first input the 
system configuration. This includes 
dividing the rotor into finite elements, 
identifying the bearing locations and 
stiffmesses, the measurement locations, 
the critical or extrapolation locations 
and the imbalance forces. The only con- 
straint in dividing the rotor into finite 
elements is that measurement and criti- 
cal locations must occur at the left end of 
an element. For good modeling accu- 
racy, the rotor must be broken into 
enough elements for adequate repre- 
sentation. This usually results in more 
computational nodes than desired 
extrapolation points. To specify which 
computational nodes correspond with 
extrapolation locations, the concept of 
Virtual measurement locations is used, 


Measurement locations can be 
declared as either real or virtual. Real 
measurement locations get their data 
directly from the data acquisition sys- 
tem. Virtual locations receive their data 
from the program's computed response 
for the node corresponding with that 
location, once the optimization is com- 
plete. The computed data at other com- 
putational nodes is discarded. 
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Typical modal probe installation on a multi-bearing machine. 
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This technique reduces the amount of 
memory required for data storage, and 
also formats the extrapolated data as if it 
were measured data, thereby allowing 
the use of the same presentation rou- 
tines for both types of data. If all the 
system Parameters — rotor configura- 
tion, suppart stiffnesses, operating 
speed, and imbalance force distribution 
— are known, the responses at the 
extrapolation locations can be com- 
puted directly. However, this is usually 
not the case. The rotor parameters 
(dimensions and materials of the shaft 
and attached disks) can be obtained 
from drawings or measurements. The 
bearing stiffnesses are a little harder to 
obtain since they depend on installation 
and operating conditions. However, rea- 
sonably good approximations may be 
obtained using perturbation or other 
techniques. The use of this program, in 
conjunction with synchronous perturba: 
tion techniques to determine bearing 
stiffnesses, is discussed by Muszynska, et 
al (1992). 

If the imbalance force distribution is 
known, the responses at the critical loca- 
tions can be computed directly without 
any optimization. However, this is rarely 
the case. The goal of the optimization 
technique is to get a good enough 
approximation to the actual imbalance 
force distribution that it can be used to 
compute the responses at the critical 
locations. Any known imbalance forces 
can be input. However, imbalance forces 
specified at this time for locations which 
will be defined as optimization locations 
will only be used to compute initial 
conditions. 

The second task shown in the flow 
diagram is to determine which elements 
contain imbalance forces that will be 
involved in the optimization process. 
Allowing control over the number of 
optimization variables lets the user trade 
length of time to get a solution with the 
resolution of the approximation to the 
imbalance distribution. At this point in 
the program, all of the initial data avail- 
able to the user has been input to the 
program. The program also uses tran- 
sient machinery data captured by a data 
acquisition system. For brevity, the opti- 
mization process will be explained at 
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only one speed, but it is actually per- 
formed at each sampled rotative speed 
contained in the data. 

The first step in the optimization loop 
is to obtain the measured data that will 
be used as the reference data. [nitial 
values are then assigned to the imbal- 
ance forces selected as optimization 
variables. The next step is to compute 
the system stiffness matrix, 

To build an algorithm that can operate 
efficiently on a personal computer, 
some limitations and special techniques 
have been implemented in the construc- 
tion of both the element and system 
stiffmess matrices. The modeling ele- 
ments must be axisymmetric cylindrical 
elements or external disks attached to 
the shaft elements with no gyroscopic 
effects, and rigid foundations at the sup- 
port locations are assumed. These lim- 
itations are not inherent to the finite 
element calculation technique. They 
have been imposed on this program to 
drastically reduce the computational 
time needed to produce the system 
response to a particular imbalance cdistri- 
bution without significantly affecting the 
accuracy of the results for the rotor sys- 
tem usec as a test case. This produces a 
much more efficient environment for 
research into optimization techniques, 


the primary research subject of this 
paper, since results from changes to the 
program algorithms can be obtained 
faster. 

Once the optimization algorithms 
have been finalized, the limitations in 
the finite element modeling can be 
removed, allowing more complex rotor 
systems to be accurately modeled. 
These limitations do make the current 
program inadequate for some 
machines. However, it is still applicable 
to a large number of rotating machines 
currently in service. This is especially 
true ifthe computed response is allowed 
the same magnitude of error usually 
inherent in the measurement process. 
These limitations force the system stiff- 
ness matrix to be the narrowest band 
symmetric formulation possible, which 
not only allows the use of the simplest 
and fastest matrix inversion algorithms 
but also reduces the amount of com- 
puter memory required to store the 
matrix. This is an important consicdera- 
tion when using personal computers 
with DOS operating systems, 

In addition to element modeling lim- 
itations, superposition is used to reduce 
computational time. Each element stiff- 
ness Matrix is constructed in layers (Fig- 
ure 3). The base layer contains all of & 











Figure 2 
Flow diagram for computer program to extrapolate measurements to other axial locations. 
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the stiffness terms for the element that 
do not contain optimization variables. If 
the element is not selected as an optim- 
ization location, the base layer contains 
all three layers, mass, shaft stiffness, and 
bearing stiffness, and the optimization 
laver is empty, 

If the shaft stiffness of the element is 
selected as an optimization parameter, 
then the shaft stiffness layer is placed in 
the optimization layer instead of the 
base laver. The same is true for the bear- 
ing stiffness layer if bearing stiffness is 
selected as an optimization parameter. 
The mass laver is alwavs placed in the 
hase layer since mass can not be selected 
as an optimization parameter. The svs- 
tem matrix is also constructed in layers, 
a base laver and an optimization layer 
(Figure 4). 

The system base layer is generated 
using the element base layers, while the 
system optimization layer is generated 
from the element optimization layers, 
This speeds Computation, since only the 
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terms associated with the optimization 
variables need to be recomputed for 
each iteration, instead of the whole svs- 
tem matrix. Note that the imbalance 
force matrix is constructed in lavers in 
the same manner as the system stiffness 
matrix. Once the force and stiffness 
matrices exist, the response matrix can 
be solved. 

Note that, for data extrapolation, all of 
the optimization variables are contained 
in the force matrix. If the system matrix 
in a solution technique can be modified 
once and then used in its modified form 
for future calculations, a significant 
reduction in computational time will 
result. This can be accomplished with 
either full inversion or LU decomposi- 
tion. This program uses a variation of LU 
decomposition, the square root method 
(Al-Khataji, et al, 1986), which takes into 
account the banded symmetric diagonal 
form of the stiffness matrix to reduce the 
computation necessary to produce the 
decomposed form. 
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Figure 3 


Layered construction of element stiffness matrix. 





System forced response equations, showing layered 
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Figure 4 


Kd= 


After the response has been calcu- 
lated, the least squared error is com- 
puted between the measured and 
calculated data. If the error is acceptably 
small, the required calculated responses 
are placed in the virtual measurement 
locations and the program proceeds to 
the next rotative speed. Ifthe error is too 
large, the program uses Powell's 
methad of discarding the direction of 
largest decrease (Press, et al, 1986) to 
control the optimization directions and 
Brent’s method (Press, et al, 1986) for 
finding the minimum along a direction 
to refine the force distribution. A new 
response matrix is computed using the 
new force distribution, and the error 
recomputed. This process is repeated 
until the error between the calculated 
and measured responses is acceptably 
small. 


Experimental results 
To test the validity of the technique 
and the program, a small rotating 
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‘construction of stiffness and force matrices. 
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machine was constructed using oilite 
bearings and a 0.375 inch (9.5 mm) 
diameter shaft 18 inches (457 mm) long 
with two heavy discs attached. The sys- 
tem was driven by a 0.1 hp (75 watt) 
motor through a flexible coupling. Syn- 
chronous measurements were taken 
with an eight-channel data acquisition, 
filtering system and then downloaded to 
the personal computer for postprocess- 
ing. Figure 5 is a block diagram of the 
system. The rotor was modeled by 
twenty-six elements (Figure 6). 

Notice that each probe location has a 
corresponding computational node 
which can be directly compared with the 
measured value during the optimization 
process. The bearing stiffness values 
were determined using synchronous 
force perturbation methodology and 
this program with the bearing stiffness 
Parameters as the optimization vari- 
ables. Once these values were deter- 
mined, they were assumed to be 
constant, which fixed all the terms in the 





svstem stiffmess matrix. At this point, a 
startup was performed and data was cal- 
lected from all of the measurement loca- 
tions. Different sets of modal probes 
were selected at each end of the rotor, 
and the responses at the other mea- 
sured locations were then extrapolated 
and compared with the actual measured 
responses. If one of the modal probes is 
located near an antinodal point, the 
extrapolated and measured responses 
are virtually identical. 

As both the modal probes approach a 
nodal point, the error in the extrapo- 
lated responses tends to increase (Fig- 
ure 7), but is still within 20 percent, an 
acceptable range for machinery diagnos- 
tic work. The noise or uncertainty of the 
measurements seems to be fairly con- 
stant along the shaft, which produces a 
good signal-to-noise ratio (SNR) near 
the antinodal points and a poorer SNR 
near the nexlal points. The quality of the 
extrapolation is directly proportional to 
the SNR. This indicates that the tech- 
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nique is highly accurate as long as good 
SNRs are maintained in the measured 
data. The computed force distribution 
required to optimize the error term is 
almost identical with the imbalance 
force used to generate the data (Figure 
8). This indicates that the method might 
enhance balancing techniques, requir: 
ing no calibration weight runs to get the 
imbalance distribution. The influence 
vectors are effectively calculated from 
the system stiffness matrix, so manual 
calibration runs can be eliminated. 


Final remarks 

The technique of extrapolating vibra- 
tian response from measured planes to 
other axial locations on the rotating 
machine, using finite element modeling, 
optimization theory and actual mea- 
sured responses, provides useful infor- 
mation for machinery protection and 
diagnostics. It works acceptably well for 
the laboratory experiments performed 
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Block diagram of experimental rotor rig and instrumentation used to 
verify operation of computer program to extrapolate vibration response. 
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Figure 6 


Experimental rotor rig shaft showing finite element construction and measurement transducer locations. 
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to evaluate whether the technique 
would work under controlled condi- 
tions. These experiments were not 
designed to be field representative. 
More work must be clone to validate the 
procedure for field use. 

However, the preliminary results are 
encouraging. IF good SNRs can be main- 
tained in the field, the technique will 
probably produce acceptable results. 
Currently, measurements tend to be 
taken at the bearings, which are usually 
near nclal points. This practice doesn't 
produce optimal measured data. How- 
ever, the greater density of measure- 
ments, when modal probes are used, 
increases the database and improves the 
technique’s performance. Additional 
data processing and different optimiza- 
tion algorithms may improve results 
even further. Both items are currently 
under investigation, Computational 
complexity can be reduced to give 
acceptable performance on personal 





computers. As computational power is 
increased, either speed or modeling 
complexity can be increased. In short, 
the procedure can enhance future diag- 
nostic systems, but further research into 
error terms and SNR increases may be 
necessary to allow the technique to func- 
tion adequately for inexperienced users 
in field applications. 


This is a subject that I am very 
interested in. 1 would appreciate 
hearing your thoughts or experi- 
ences on this subject. I would be 
grateful for any information you 


can send me, including actual 
machine data and case histories. I 
will keep you informed of the pro- 
gress of my work. 


Donald FE. Bently 
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Extrapolated data for measurement planes 3, 4, 5, and 6 using measurement planes 1, 2, 7, and 8 as modal measurement locations, 
presented in the form of Polar plots of synchronous vibration. Solid lines are measured data, broken lines are extrapolated data. 
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Figure 8 





Imbalance force distribution versus rotative speed, created during optimization process. Solid lines are computed responses. Broken lines 
are the imbalance forces used to create the transient data. The left plot is the direct component of the imbalance force, while the right is the 
quadrature component. “Direct” and “quadrature” refer here to orthogonal coordinates referenced to the measurement transducers. 
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